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Abstract 

We consider a simple model for the coordinate-space vacuum polarization function which is often 
parametrized in terms of a screening mass. We discuss the circumstances in which the standard 
result for the screening mass , m sc = ttT, is obtained. In the model considered here, that result 
is obtained when the momenta in the relevant vacuum polarization integral are small with respect 
to the first Matsubara frequency. 
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In a number of recent works [1-3] we have calculated various hadronic correlation func- 
tions and compared our results to results obtained in lattice simulations of QCD [4-6]. The 
lattice results for the correlators, G(r, T), may be used to obtain the corresponding spectral 
functions, a(cu,T), by making use of the relation 

PCX) 

G(t,T)= dcua(cu,T)K(r,u,T) , (1) 
Jo 

where 

K , , , T , cosher - 1/2T)] 
sinh(u;/2T) 

The procedure to obtain a(cu, T) from the knowledge of G(t, T) makes use of the maximum 
entropy method (MEM) [7-9], since G(r,T) is only known at a limited number of points. 

In our studies of meson spectra at T = and at T < T c we have made use of the Nambu- 
Jona-Lasinio (NJL) model. The Lagrangian of the generalized NJL model we have used in 
our studies is 



C = q{ij - m°)q + ^ ^[(gA'g) 2 + [qil^qf] (3) 

i=0 

" ^E[(^7^) 2 + (^ 7 57^) 2 ] 

i=0 

G 

+ -^{det[g(l + X 5 )q] + det[g(l - A 5 )g]} + C conf . 

Here, m° is a current quark mass matrix, m° = diag(m Q u) m Q d) m° s ). The \ are the Gell- 
Mann (flavor) matrices and A = a/2/31, with 1 being the unit matrix. The fourth term is 
the 't Hooft interaction and £ con / represents the model of confinement used in our studies 
of meson properties. 

In the study of hadronic current correlators it is important to use a model which respects 
chiral symmetry, when m° = 0. Therefore, we make use of the Lagrangian of Eq. (3), while 
neglecting the 't Hooft interaction and C con f. In order to make contact with the results of 
lattice simulations we use the model with the number of flavors, Nf = 1. Therefore, the A* 
matrices in Eq. (3) may be replaced by unity. We then have used 

C = q{ij - m°)g + ^-[{qqf + (to) 2 ] (4) 

- ^r[(?7M<?) 2 + (9757m?) 2 ] 
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in order to calculate the hadronic current correlation functions in earlier work [1-3]. 

In order to present our results in the simplest form, we consider only the scalar interaction 
proportional to (qq) 2 . We also extend the definition of <j(uj,T) of Eq. (1) to include a 
dependence upon the total moment of the quark and antiquark appearing in the polarization 
integral. Thus we consider the imaginary part of the correlator, a(oj, P). Since we place P 
along the z-axis this quantity may be written as cr(u, 0, 0, P z ) in accord with the notation 
of Ref. [10]. In this work we will present our result for the coordinate-dependent correlator 
C(z) which is proportional to the correlator defined in Eq. (1) of Ref. [10], 

a(u,0,0,P z ) 



C{z) 

We may also use the form 

C(z) = 



dP z e lP * 



du- 
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We have made a study of the screening mass in a simple model in order to understand 
the origin of exponential behavior for the correlator. To that end we make use of Ref. [11]. 
We consider the Matsubara formalism and note that the quark propagator may be written, 
with = 1/T, 



Sp( k,u n ) 



7°(2n + l)?r//3 + ~j ■ k - M 



(2n + 1) 2 tt 2 //? 2 + k 2 + M 2 
For bosons the vacuum polarization function is given as Eq. (1.51) of Ref. [11], 

d 3 k 1 1 



(7) 



n (^°) = §E (27r) 3 4nV 
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(8) 



+ p°) + ( k + ~p ) 2 + M 2 



We modify Eq. (8) to refer to fermions. In this case the Matsubara frequencies are 

(2n+ 1)tt 



UJr, 



(9) 



and we have 
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d 3 k 



(-f°7l/P + ^ ■ ~k ) ( 


7°(p° + 7r//3) + 7^ • (k + 
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FIG. 1: The function C(z) of Eq. (6) is shown for a sharp cutoff of k max = 0.1 GeV. The dotted 
line represents an exponential fit to the curve using m sc = 1.23 GeV. (We recall that irT is equal 
to 1.27 GeV.) 

if we keep only the first term in the sum, where u = ir/fl. As a next step we drop p°, so 
that we have 



n(?,o) = f> 1 <n ' 



(11) 



2(3 J (2tt) 

We then take p along the z axis and write H(p z ) = n( p , 0). We define 

C(z) = J dp z e^ z U( Pz ). (12) 

In our calculation we replace g 2 /2/3 by unity and use a sharp cutoff so that | k \ < k max . 

The results of our calculation of C(z) of Eq. (6) are given in Figs. 1 and 2. In Fig. 1 we 
use k max = 0.1 GeV and in Fig. 2 we put k max = 0.4 GeV. For our calculations, we have 
m sc = tiT = 1.27 GeV when T = 1.5 T c and T c = 0.27 GeV. Thus, the k max values considered 
here are less than m sc and that feature leads to the exponential behavior seen in Figs. 1 
and 2. If k max is made larger than 0.4 GeV we begin to see deviations from exponential 
behavior for C(z). (Since in our calculations reported in Refs. [1-3], the integrals were 
regulated with a Gaussian regulator exp[— k 2 /a 2 ] with a ~ 4 GeV, we can see that the k 
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FIG. 2: The function C(z) of Eq. (6) is shown for a sharp cutoff of k max = 0.4 GeV. The dotted 
line represents an exponential fit to the curve using m sc = 0.961 GeV. (We recall that irT is equal 
to 1.27 GeV.) 

values in those calculations are so large as to preclude obtaining exponential behavior for 
our coordinate-space correlator.) 

The goal of this work was to consider a simple quark model for the calculation of a 
hadronic current correlation function and to determine the conditions under which the 
coordinate-space correlator is dominated by the screening mass which is given by the first 
Matsubara frequency. We have found that the standard result is obtained if the quark 
and antiquark momenta in the vacuum polarization calculation are small compared to that 
frequency. 
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